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We suggest a theoretical description of the force-induced translocation dynamics of a polymer 
chain through a nanopore. Our consideration is based on the tensile (Pincus) blob picture of a 
pulled chain and the notion of propagating front of tensile force along the chain backbone, suggested 
recently by T. Sakaue. The driving force is associated with a chemical potential gradient that acts on 
^SJ , each chain segment inside the pore. Depending on its strength, different regimes of polymer motion 

(named after the typical chain conformation, "trumpet", "stem-trumpet", etc.) occur. Assuming 
c 2 ^ ■ that the local driving and drag forces are equal (i.e., in a quasi-static approximation), we derive an 

04 \ equation of motion for the tensile front position X{t). We show that the scaling law for the average 

■ translocation time (r) changes from (r) ~ ]<l'^^ j f'^l^ to (r) ~ N^'^" / f (for the free-draining case) 

D \ as the dimensionless force fR = aN" f /T (where a , N , v, f, T are the Kuhn segment length, the 

pL^ ■ chain length, the Flory exponent, the driving force, and the temperature, respectively) increases. 

These and other predictions are tested by Molecular Dynamics (MD) simulation. Data from our 
computer experiment indicates indeed that the translocation scaling exponent a grows with the 
pulling force fn) albeit the observed exponent a stays systematically smaller than the theoretically 
predicted value. This might be associated with fluctuations which are neglected in the quasi-static 
i ' approximation. 
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I. INTRODUCTION 



& 

(~| ' The force-induced translocation through a pore in the membrane is mainly motivated by the possibility for fast 
O ■ DNA and RNA sequencing P, |3| ■ Translocation dynamics is an essential component of transport in biological cells 
^ . Most of the experimental studies so far deal with a driven polymer translocation which is realized by applying an 
electrical field across a narrow pore JJ. However, despite extensive research, numerous computer experiments, and 
^sj ■ a variety of attempts for theoretical interpretations , the translocation dynamics even in the generic case remains at 
J> ] present insufficiently well understood [5l-[l2|. As emphasized in a recent review paper [l^, currently there exists a 
■ plethora of theoretical predictions for the value of the translocation exponent a which governs the scaling of the mean 
translocation time (r) with the length N of the polymer chain ((r) cx iV"). These values vary in a rather broad interval: 
one predicts a = 1 [3, H, [1], or a = -I- 1 — 7i w 1.49 jlO| for a three-dimensional {3D) self-avoiding (SAW) polymer 
chain (where the Flory exponent s» 0.59 and the surface exponent 71 w 0.68), up to a = (1 + 2i^)/(l + ^ 1.37 
O; [ni anda = l-Fi/R:!l.59 [i[llj. 

In a recent paper we demonstrated that the dynamics of a single polymer moving through a nanopore in a 
membrane (i.e. the translocation dynamics) with or without an external force being present can be treated within 
• • ■ the framework of fractional Brownian motion (fBm). It was shown that the corresponding Fokker-Planck equation 
. ^ [ of motion (FPEM) contains time-dependent drift and diffusion terms. In the case of non-driven chain this FPEM 
^> . leads naturally to anomalous diffusion of the translocation coordinate s{t) [13] that describes the number of segments 
^ ' threaded through the pore at time t. On the other hand, for a driven translocation fluctuations play a relatively 
, moderate role and one can directly construct reasonable scaling relations that describe the process. Recently, an 
interesting approach based on the notion of tensile force pro pag ation along the chain (within the well known picture 
of tensile blobs [20I [2l| ) has been suggested by T. Sakaue |l5l. [l6j in an effort to provide a consistent description 
of translocation dynamics. Sakaue's approach appears physically sound and well capable of providing a plausible 
interpretation of existing observation. However, his mathematical treatment requires the use of a "cut-off" trick as 
applied to the polymer segments density function (see more details in the Appendix . This trick is questionable 
and affects the predicted translocation exponent a. 

In the present paper we follow Sakaue's approach and consider theoretically the case of a driven translocation by 
means of a different derivation which does not require additional conjectures. In Sec. II we investigate the various 
regimes (named traditionally a "trumpet" , "stem-trumpet" , and "stem" regimes) , observed for weak, intermediate 
and strong driving forces. Our consideration is based on equating the driving and drag forces, i.e., on a quasi-static 
approximation whereby fluctuations in this balance of forces are ignored. Our main concern is the scaling law for the 
mean translocation time (r) with respect to the chain length N and the applied driving force /. We find (t) oc N°'f~^ 
albeit with a different expression for a as well as the mean translocated length M{t) vs. time. 



2 



We also perform extensive Molecular Dynamics (MD) simulations in order to test our theoretical predictions. In 
Sec. Ill we describe the simulation model and compare theoretical predictions with simulation data. This comparison 
shows that the exponent a, found in MD-simulation, is systematically smaller then its theoretical estimate even 
though a changes with / as expected. We argue that this may be due to fluctuations which are neglected in the 
theoretical treatment. Finally, in the Appendix |^ we explain in detail why the calculations of Sakaue (l5l . [l6j are 
bound to produce a scaling exponent that differs from the proper one. 

II. DYNAMICAL SCALING PICTURE OF A DRIVEN TRANSLOCATION 

A. A general overview 

We assume that the pulling force is applied only to a bead which is inside the pore, i.e., the potential drops mainly 
across the pore. This assumption was recently examined for charged DNA translocations J-7, 18]. It was shown, in 
particular, that the ionic current generates a nonuniform electric field which is extended far beyond the pore and acts 
on both the DNA and its counter ions. In this case the DNA threading through the pore is preceded by its capture by 
electroosmotic flow. On the other hand, it can be shown that for a sufficiently narrow pore the voltage drop mainly 
takes place inside the pore ^1^. In this situation the translocation is mainly determined by the threading process 
which is the main subject of our investigation. 

When a pulling force is suddenly switched on, tension starts to propagate along the chain backbone which alters 
the polymer conformation progressively. Eventually, after some characteristic time a steady state is reached and the 
whole polymer starts moving with constant velocity. Such a non-equilibrium response is of great importance in many 
biological and technological situations. Below we follow to some extent the model developed in the recent papers by 
Sakaue 0,13. 




FIG. 1: (Color online) Dynamical response of the polymer chain shortly after the driving force / has been switched on. By 
the time t the tension is transmitted up to the N{t)-ih monomer. The number of the translocated monomers is denoted by 
M{t) (indicated by a dashed line on the trans-side of the pore). The distance between the propagating tension front and the 
membrane is marked as X{t). The portion of the chain marked as X{t) is moving with an average time-dependent velocity 
v{t). a) The "trumpet" regime at driving force kBT/aN'^ ^ / ^ ksT/a . The initial blob size is ^in = ksT/f. b) At 
ksT/a ^ f <^ {kBT/a)N'^ the part of the chain affected by tension starts as a "stem" of the length S{t) and then turns into 
a "trumpet" ("stem-trumpet" regime), c) At / > {kBT / ojN" the tensed part of the chain is completely stretched ("stem" 
regime). R marks the unperturbed chain size, i.e., R — aN'^ . 

Figure [1] shows schematically the case of tensile force transmission. The driving as well as the friction forces are x- 
dependent which leads for a relatively small driving force exerted on a segment in the pore, ksT/aN'^ ■C / ^ ksT/a, 
to the so-called "trumpet" blob picture (shaded area in Fig. [T^) jl^, [2lj . By equating the local driving and drag 
forces (see Eq. (|2.5p ) and using the relationship between the tensile blob size ^(x) and the force f{x) (cf. Eq. (|2.4p ). 
one could derive the "trumpet" profile equation (cf. Eq. (|2.9p ). The location of the tension front X{t) is defined 
by the free boundary condition f{x = —X{t)) = 0. This condition along with the material balance Eq. (|2.14l) and 
the closure relation, Eq. (|2.15l) . leads to the differential equation, Eq. (|2.20p . for the front propagation X{t). The 
characteristic time of this transmission is given as ri = Ci iYi+^Z/^-i-i/i^ (cf. Eq. (12231)), where the dynamical 
exponent z = 2 -\- 1/ v for the free-draining (or Rouse) case, and z — 3 for the non-draining (or Zimm) case, Ci being 
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a model dependent constant. For a Rouse chain, ri = Ci N^'^^ / f which agrees with the scaUng law predicted by 
Kantor and Kardar [l^]- During the process of tensile force transmission (or, front propagation), the velocity v{t) 
of the moving domain decreases (because as time goes by more and more segments get incorporated in the moving 
domain) and at the time t = ti it approaches a stationary value Vs- After that a stationary regime sets in and the 
rest of the chain is sucked into the pore with a constant velocity Vs (cf. Eq. (|2.24p ). The characteristic time of this 
stationary suction process is T2 — C2A^^'^//^~^, (see Eq. ()2.25p ). The times n and T2 add up to a net translocation 
time (r) = ti +r2. Generally, the force transmission time ri prevails for a strong enough driving force and a long chain 
(see Eq. (12.271) ). Thus, the translocation time scaling ranges between iV^'^//^~^ and iV^+''//^~^~'^/'^, depending on 
the pulling force intensity. 

For moderate, ksT/a ^ / ^ {ksT /a)N'^ , and strong, [ksT /a)N'^ < /, forces, the "trumpet" regime is replaced 
by the "stem-trumpet" , see Fig. [T]d, and "stem" , see Fig. [TJ;, regimes respectively. In the "stem-trumpet" case, the 
part marked as S{t) of the total moving domain looks like a stem whereas the rest resembles a trumpet shape. For 
the "stem" scenario of the force transmission the total moving domain looks like a completely stretched portion, see 
Fig. [TJ;. In both regimes the characteristic times of the force transmission and stationary suction into the pore are 
n ~ N^^'^ /f, and T2 ^ N'^'^/ f. It could be seen that the time ti dominates for a longer chain (see Eq. (|2.45p ). In 
the TableUwe summarize the results of the theoretical scaling prediction for the translocation time (r) oc N^f^^ (for 
the Rouse case) as well as the corresponding MD-simulation findings which will be discussed in this Sec. [ll] and in 
Sec. mi 



Regime 


Exponents a and S (theory) 


Conditions 


Exponents a and S (MD) 


Trumpet: 


a = 2;/ ^1.18, 5 = 1.66 
a = l + u^ 1.59, 5 = 1 
sec Ea. (|2.26|) 


Ir < (Ca/Ci)"'''^-''' 
fR > (Ci/CiY^^'-'^ 
see Eq. (12.271 


a ^1.11, (5 « 1.17 
a « 1.47 , S ^ 0.97 
see Sec. [m] 


Stem- Trumpet 
and Stem: 


a = 2uP^ 1.18 , (5 = 1 
a^l + v!^ 1.59 ,5 = 1 
see Eq. (|2.44|) 


short chain 
longer chain 
see Eq.(|2.45|l 





TABLE I: Exponents a and S for the translocation time scaling (r) oc TV"/ in the free-draining (Rouse) case as predicted by 

the theory (theory) and by our MD-simulation (MD) (see Sec. The parameter /r =* aN" f /ksT , whereas Ci and C2 are 
some numerical model-dependent parameters. 

Schematically, our theoretical findings for the scaling behavior of (r) for the "trumpet" regime (see Sec. IIC) are 
plotted in Fig. [2^ (for (r) ~ N"" ) and in Figl^b (for (r) ~ f^^)- The time evolution of the number of translocated 
monomers M{t) is also investigated in Sec. IIC2. Again, depending on driving force intensity, the M{t) vs. time 
t scaling law changes, see Figl2}:, from M{t) oc (/t)^/(^+'') (at strong force, when the force transmission dominates) 
(see Eq. p.30|) ) to M{t) oc {ft)/N'^ (for weak force, when the stationary suction process dominates and t ^ T2) (see 
Eq. (I2.34p ). These M{t) vs. t dependencies clearly illustrate the presence of an initial stage, which in the Sec. IIB is 

called blob initiation, and is confirmed by our MD-results (see Fig. [7] where the translocation coordinate (s) M{t) 
vs. time t is given). It will also be shown that the characteristic time of this initial stage Tjnit changes with force as 
Tinit 1// (see Eq. ^ and Fig. E 



<t:> 




FIG. 2: Schematic representation of the main scaling predictions: a) Translocation time (r) vs. chain length A'^ at a fixed force 

/. With growing force fR =' aN" f /ksT the scaling law (r) ~ N°' changes from (r) ~ N'^'^ to (r) ~ A^^". b) Translocation 
time (r) vs. driving force / at different chain lengths Ai > A2. c) Translocation length M{t) vs. time t for two different forces 
/2 > /i. The characteristic time for initial blob formation goes as rinit ~ 1//. 

We will devote the rest of this Sec. [ll] to the detailed discussion of the "trumpet" , "stem-trumpet" and "stem" 
regimes as well as to the mean translocation length M{t) behavior. Readers who are not interested in the mathematical 
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aspects of the problem can skip these details and go to Sec. IIIII One should keep in mind, however, that before one 
or another type of tensile force transmission comes into play, an initial blob should always set in. 



B. Blob initiation 



The starting deformation (i.e., blob initiation) of a polymer chain that has been initially at rest and is then pulled 
by a force / constitutes a necessary stage to set in the driving force transmission. An initial blob of size ^init = ksT/ f 
is formed which contains g beads so that ^init = ag'^ where a is the Kuhn segment length and v is the Flory exponent, 
mentioned above. Denote tq = a'^C,o/kBT as the characteristic diffusion time of a single bead with Co being the Stokes 
friction coefficient [22i] . Then the blob initiation time can be written in the scaling form 

r,,t = ro0(^^) (2.1) 

so that the scaling function (/>(j/) depends only on the dimensionless combination ^mit/a (i-e., on the size of the blob 
itself). On the other hand, we assume that the initial blob formation happens due to purely mechanical traction which 
almost docs not affect the chain meanders. This is supported by the fact that within the time interval < t < rinit 
only few segments are translocated (see Fig. [2t). Thus, the characteristic time Tjnit does not depend on temperature 
T and is determined only by the Stokes friction coefficient Co and the pulling force /. Since tq oc l/ksT, this implies 
(f>{y) oc y, and one gets as a result 

a\o\ fkBT\ Coa „x 
oc — — l^-^) 



ksTj \ af J f 

The rearrangement of the initial blob leads a.t t > rinit to a subsequent tensile force transmission which is governed by 
the local balance of driving and drag forces. The blob initiation characteristic time rinit is indeed clearly seen in our 
MD-simulation results (see Figl?]) and found to change as ^ 1//. In the next sections we consider different regimes 
of this tensile force transmission, depending of the strength of /, which are named after some characteristic shapes of 
the resulting chain conformation. 



C. Trumpet regime 



In this case the driving force is moderately strong and falls within the range 

ksT ^ ksT 



(2.3) 



Schematically, a typical conformation of a driven polymer in the trumpet regime is represented by the shaded area 
in Fig. [TJ At a given time t some fraction of the chain is subjected to tension. This part of the chain comprises a 
sequence of Pincus blobs of different size. 



e(x) 



keT 



(2.4) 



where f{x) is the local (i.e. x-dependent) value of the driving force. As is evident from Fig. [U the blob size at a 
given position x corresponds to the lateral chain meandering at x. According to the blob definition ^(x) — ag{x)'^ 
where g{x) is the number of beads in the blob (we recall that the chain inside the blob follows Self Avoiding Walk - 
statistics). On the other hand, as stipulated by the local force balance, the driving force equals the friction (Stokes) 
force, i.e. 



fix) = Co vit) 



X 

1 






a 



-X(t) 



dx' 

W) 



(2.5) 



In Eq. (|2.5|) we take into account that dx'/S,{x') is the number of blobs on the interval {x',x' + dx'} whereas 
Co v{t) [^{x')/aY^'^ is the local Stokes friction force. The dynamical exponent z in Eq. ()2.5p . which determines 
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the total number of monomers in this portion of the chain, is z = 2 + 1/v for the case of Rouse dynamics, and z = 3 
for the Zimm statistics [2l|. Thus, taking into account Eq. (I2.4p . the expression for the blob size becomes 



or, represented in a differential form. 



Co v{t) 



Eq. ()2.7p should be supplemented by the boundary condition 



/ 



(2.6) 



(2.7) 



(2.8) 



which represents the initial blob in the 'trumpet' (cf. Fig. [T]). Then, the solution of Eq. (|2.7p takes on the form 

-21 1/(2-^) 



^(x) = a 



Co ^(0 
kRT 



kgT 
f 



knT 



(2.9) 



where in the case of Rouse dynamics 2 — 2 = 1/;^. 

One can fix the blob size by the requirement that the force at the free boundary, i.e., at a; = —X{t), is equal to 
zero. Therefore, 



fix = ^X{t)) 



knT 



0, 



^x = -X{t)) 

or, in view of Eq. (HH), this yields (Co v/kBT)(~X{t) + keT/f) + {af/kBTy~^ = 0. Solving the latter for X{t), 



(2.10) 



X{t) = 



ksT ksT f af 



f Cov{t) ykeT 
one may exclude / in Eq. (|2.9p in favor of X{t). As a result we have 



z-2 



[Covit) ix + X{t))/kBTf^'-^^ 



(2.11) 



(2.12) 



The largest blob of size ^x{t) in the trumpet is placed a.t x = —X{t) + ^x{t)- From Eq. (|2.12p one immediately 
obtains 



knT 



Co a v{t) 



(2.13) 



We are now in a position to derive the differential equation for the location of the tension front X(t). To this end we 
use the material balance equation (mass conservation law). Let the tension be transmitted up to the N(t)-t]i monomer 
whereas M{t) denotes the number of translocated monomers (see Fig. [T]). Then the material balance equation reads 



1/^ 



-X(t) 



dx 

W) 



M{t) = N{t) 



(2.14) 



where the first integral yields the number of monomers in the shaded portion of the chain (see FigH]). 

Denote the total number of monomers that has been subjected to tension during the time interval {0,t} by N{t). 
At t = all these N{t) monomers have been in equilibrium occupying a region of average size X{t). Therefore N{t) 
and X{t) are related by the Flory expression, i.e.. 



X{t) = aN{tY, 



(2.15) 



which will be used below as a closure relation. 
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The local density in the moving domain of blobs (shaded portion in Fig. [T]) is p{x) ~ 9{x)/^{x)^ whereas the 
cross-section area is locally ~ The smallest blob next to the membrane orifice seeps through the pore 

into the trans-side of space so that the number of translocated monomers can be calculated[l^ as 



Mit) = J p{xMx)\^^,v{t')dt' = 



t 

= - / mf-^ \xit') 

To J L . 



ax = 0) 



dt' 



i[x = 0) 



w(t') dt' 



(2.16) 



where we have used Eq. (|2.12p for the blob size and have introduced the bead characteristic time tq = a^C,Q/kBT as 
well as the dimensionless (tilded) variables (TB]: v{t) = (^oav{t)/kBT and X{t) = X{t)/a. Moreover, using Eq. (j2.11|) 
we can exclude v{t') in favor of X{t) and the dimensionless force fa — af/kBT. As a result 



To J 



dt' 



Xif) 1-1/faXit') 



(2.17) 



Now we use Eqs. Eq. (|2.11l) and (I2.12p to determine the first integral-term in Eq. (|2.14p where we also take into 
account Eq. ()2.17p and the closure Eq. (|2.15p . As a result, the material balance condition Eq. (|2.14|) takes on the 
form 



X{t) 



1 - 



1 



faXit) 



~l-L ~z-l-i 1 
fa " fa " 

To 



dt' 



X{t') 



faX(t') 



(2.18) 



For relatively long time intervals faX{t) fa — 1 and Eq. f|2.18p takes on the form 



x{t) fi + /: ■■ 



dt' 



X{t) 



To J Xit') 



Differentiation of Eq. (|2.19p with respect to time yields a differential equation for X{t) 



(2.19) 



1 - Bo faX{t) 



dt X{t) 

where Bq is a constant of the order of unity. Eq. (j2.20p should be supplemented by the initial condition 

X{t = 0) = Cin = i 
J a 



(2.20) 



(2.21) 



which simply indicates that the force transmission starts right after the initial blob has formed as we have discussed 
in Sec. II A. The solution of Eq. ()2.20|) has the following form 



t = to + ToBofl'''-'-+' X{tY/''+^{l-Col faX{t) } 



(2.22) 



where Co = I/-B0 and to = tq{1 — Bo)/ fa- The characteristic time ri for the last monomer to attain its steady state 
follows from X{ti) ^ N" which, due to Eq. (|2:22|) . yields 



n - TO Bo /y--+i^i+- [1 _ Co/rj"-' 

where N is the total chain length, fji = Rf/kBT and R = aN^ stands for the unperturbed chain size. 



(2.23) 
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1. Stationary part of the translocation: suction into the pore 

The stationary regime sets in within a characteristic time ti. After that the non-translocated part of the chain, 
N — M{ti), is puhed as a whole entity towards the pore (suction). The stationary velocity Vg is defined from Eq. 
((2ll|) by putting X{ti) = N" . As a result for JaX{Ti) > 1 we have 

^s-^ J:-' (2.24) 

where the dimensionless velocity Vs = Co ciVs/kBT. The characteristic time for the stationary part of translocation is 
then derived as 



T2 = TO -4r^ = ^0 (2-25) 



It 

Eventually, the total translocation time (t) can be seen as a sum of ri and T2. Thus, by taking into consideration 
Eq. (|2.23|) (at fa ^ 1) and Eq. (|2.25|) we derive one of our main results, 

(t) = TO CJl'''-'+^N^+'' + TO C2 Jl-' N^" (2.26) 

where Ci and C2 are some numerical model dependent constants. The first term in Eq. (|2.26p dominates in case 

C'l/a^""^^^-^^^" ^ C2 /a~^ N'^'^ , that is, for sufficiently large driving forces / and moderate chain lengths or, 
alternatively for very long chains even if the driving force is rather weak: 

Jr » (C2/Ci)''/(i-'') (2.27) 

where fn = aN^ f /ksT . The combined schematic representation of the scaling law, Eq. (|2.26p . as well as the criterion 
given by Eq. (I2.27P is given in Table |T] and in Fig. [2^, b. 

The central result given by Eq. (|2.26l) predicts a scaling relation for the mean translocation time (r) oc N°'f~^ with 
a force-dependent exponent a that grows from 2:/ to 1 -I- as the parameter aN'^ f /ksT increases. Numerical results 
obtained by means of Langevin Dynamics (LD) simulation in two dimensions (2D) clearly support this behavior: one 
observes a = 1.50 for relatively short chains {N ^ 200) and a = 1.69 for longer chains I23'|. In three dimensions (3D) 
there is no clear evidence of such crossover. In the total interval of chain lengths which was studied in 3D the exponent 
a — 1.41 [2^ and a = 1.36 [l^l , i.e. the exponent value falls in the range between 2j/ and 1 + i^. In a series of LD 
simulation studies by Lehtola et al. (25l - [27j it was found (within chain length interval N ^ 800) that a systematically 
increases with / up to roughly a — 1 + v. This is in agreement with our theoretical and MD-results (see below Sec. 
imp but is in clear contradiction to the results of K. Luo et al. [28| where the exponent Q; = l + t^«1.59 was found for 
relatively short chains, N ^ 200 and for small driving forces. Moreover, the exponent a becomes smaller, a « 1.37 
as the driving force grows. Experimentally, the scaling law (t) oc N'^'^ / f , which follows from Eq. (|2.26l) for relatively 
weak forces fn and z = 3 (Zimm dynamics), has been found in ref. [29l. [sfl in the case of synthetic pore. Recently 
the scaling given by the first term in Eq. (I2.26P has been obtained in ref. I9| using the so-called "iso-flux trumpet" 
model (which was also inspired by the Sacaue's paper). 

Eventually, we should like to point out that the first term of Eq. (|2.26p (which corresponds to the characteristic 
time of tensile force transmission) differs from the corresponding term given by Eq. (10) in ref. [16]. This difference 
arises even if one assumes, following Sakaue's idea, that the entire trumpet is characterized by a single time dependent 
velocity ( "iso- velocity trumpet"!). The reason for this discrepancy is discussed in the Appendix \K\ 

2. Time evolution of the translocated portion M{t) 

The number of translocated monomers M(t), given by Eq. (|2.17h . grows with elapsed time. In the case when the 
second term in Eq. (|2.22p dominates, one gets 

i ~ TO So II'''^'^' X(t)i/''+i (2.28) 
Combining this equation with Eq. (|2.17p (where again we assume faX{t) ^ 1), one arrives at 

t/T-O 

M{t) = {lr I (t)'^ dt^Coifar{ty/^'+''^ (2.29) 
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where the exponent x = (-^^^ — 1 — v)/v{l + v). In the case of Rouse dynamics zv = 2v ^ 1 and x = (1 + i/) . The 
resulting translocation relationship can be written as 



M{t)=Co{fat 



1/(1 + '') 



(2.30) 



As far as mean translocation time (r) should satisfy the requirement M{(t)) — N, one recovers the result (r) ~ 
roN^'^'^/f (see Eq. STm ). 

In the case of small driving force, /a^V" — 1 and faN"^ ^ {C2/Ciy^'^^^'^\ the translocation process is mainly 
determined by the stationary suction into the pore whereby the translocation time (r) ~ C2TqN'^'^ / f^~'^ . On the 
other hand, the number of translocated monomers according to Eq. (|2.16p reads 



M{t) = ' 



C(x = 0) 



1 



e(x = 0) 



dt = [i 



dt' 



TO J [Xit')]^ 



(2.31) 



where we have used the expression for the blob size Eq. (|2.12p with the stationary velocity Vg- In the stationary 
regime the size of the moving domain has the form X{t) = X{ti) — Vst/ro TV — Vst/rQ. By making use of this in 
Eq. (|2.3ip and after integration, we arrive at the result 



M{t) 



t 

T2 



(2.32) 



Consider first the large time limit t ~ r2. In this case M{t2) ~ N and the Eq. ()2.32p leads us back to the condition 
faN" ~ 1. At f ^ T2 the expansion in Eq. p.32p suggests that 



M{t) ~ ci 



fa 



'z-l-l/u 



t 



where, as before, t — t/rQ. In the case of Rouse dynamics we finally obtain 

M{t) ^ ^ /J 



(2.33) 



(2.34) 



As can be seen from Eq. (|2.34p . the formal extrapolation the linear time dependence law up to the final time tf 
(where M{tf ) ~ N) gives </ oc N^^'^ / fa ~ ti ^ T2. As a consequence, a sub-linear slowing down behavior should be 
seen at a later time interval. The two regimes mentioned above (i.e., for strong and weak driving forces) are illustrated 
in Fig. dt. 



D. "Stem" and "stem-trumpet" scenario 



1. Stem 



The above consideration is valid in the range of forces given by Eq. (|2.3p . For stronger forces which are within the 
range 

1 ^ 7a « iV" (2.35) 

the tensile force transmission follows the so called stem-trumpet scenario 0, [2l|, as indicated in Fig. [T|d, whereby 
the part of the polymer chain that is close to the membrane attains a completely stretched conformation. 

Consider first the limit of very strong force fa > where the configuration is depicted by the stem picture shown 
in Fig. [TJ;. In this case the force balance reads 

v{t)X{t) = fa (2.36) 

The material balance (cf. Eq. (|2.14p ) yields 



X{t) + M{t) = N{t) (2.37) 
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The number of translocated segments (cf. Ea. (l2.16|) ') is 



M{t) = - / v{t')dt' = — / v{t')dt' 
a J To J 



(2.38) 



With the help of Eqs. ((OS)) . (fOT)) and as well as using the closure Eq. (P?T5|) one finds 



Xit) 

In differential form this equation reads 

TO I Co 



fa f dt' 



TO J X{t') 



X{t) 



111, 



(2.39) 



X{t) - C\X{t) 



dX 
~dt 



fa 



(2.40) 



with an initial condition X{t = O) = I and Co, and Ci being some constants. The corresponding solution of Eq. 
([2^ is given by 



—t — Ci + Co — Co 
To 



X{t) 



1 - 



Ci 



Co 



(2.41) 



Again the characteristic time t\ is defined as X(t\) = R = iV". Then, for the reasonably long chains X^ " ^ Ci/Cq 
the first terms in Ea. (|2.4ip prevails. As a result, the characteristic time t\ is obtained as 



T\ = To 



fa 



(2.42) 



As before, the second (stationary) stage of the translocation is related to the pulling force by considering this part 
of the chain as an entity which moves with stationary velocity Vs — fa/N^ . The corresponding characteristic time T2 
is given by 



R N^" 

T2 = To — = To 

fa 



Therefore, the total translocation time (r) is a combination of ri and T2, i.e. 



(r) = Ci To— ^ + C2 TO 

fa fa 



The first term in Eq. p.44p dominates at 



(2.43) 



(2.44) 



(2.45) 



i.e., for a reasonably long chains. One should note that in this limit we recover the scaling law of the trumpet regime, 
provided Eq. (12:45]) holds and X{t) = [JatY ' ^^^''^ The velocity v{t) = Ja/X{t) = Ja{Jat)-''l'^^+''\ and the number 
of translocated monomers is obtained as 



M(<) = i jv{t')dt'^Bo (^fat) 



1/(1+1') 



(2.46) 



whereby one goes back to Eq. (|2.30p . 
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2. The Stem-trumpet regime 

In the interval of force strengths, Eq. (|2.35p . the chain deformation starts with the formation of a stem whereby 
the velocity changes as 



1/(1 + '') ^To\''/(l+'^) 

T 



(2.47) 



Thus, the velocity decreases with time and at some moment t ~ the drag force at the stem-trumpet junction 
becomes comparable to kgT/a, i.e., Co^^('''^) — ksT/a. Therefore, 



To [fa 



(2.4 



At t > the stem-trumpet picture dominates. For the blob size of the flower (trumpet) part the same differential 
equation Eq. p.7|) holds. However, the boundary conditions (BC) are different. Namely, in this case ^{x = —S{t)) = a 
and the solution of Eq.^J^ yields ^(x) = a{l + Cov{t)[x + S{t)]/T}^^'''^~''\ Therefore, the solution is given in the 
total cc-interval by 



a{l + Cov{t)[x + S{t)]/kBT}^^^^-'^ , for - Xit) < x < ~-S{t) 
a , for - S{t) <x<0 



(2.49) 



For X = —X{t) the force vanishes, i.e., f{x — —X{t)) = T/^{x = —X{t)) = 0, and one has 



X{t) = S(t) 



The material balance may be written as 

-S(t) 



/ 






a 



-X{t) 



1/u 



dx 1 



+ - Sit) + M{t) = Nit) 



By making use of Eqs. ()2.15p . (|2.49p and (|2.50p in Eq. (|2.5ip . one arrives at 



(2.50) 



(2.51) 



X(t) + - / v(t')dt' = X{t) 
a ' 



(2.52) 



In order to exclude the velocity from Eq. (I2.52p we should note that the tensile force a.t x — —S{t) (junction point 
between "stem" and "trumpet" ) is given by 



fix = -Sit)) 



(2.53) 



Then the force balance for the stem can be written as 

■Sit) 



Co«(t) 



and consequently. 



f-fix^-Sit))^f- 



vit) 



kuT 



(2.54) 



(2.55) 



This equation together with Eq. ()2.50p yields 



Xit) = ^^fa 

vit) 



(2.56) 
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and one retrieves the corresponding equation Ea. (j2.36p for the stem case, 
arrives at the differential equation for X{t) 



Now, using Eq. ((2?56l) in Eq. (|232)) one 



To Bo 



Xit) 



- BiX{t) 



dX{t) 
dt 



= fa 



(2.57) 



which is again equivalent to Eq. (I2.40p for the "stem" case. As a consequence, the expressions for the translocation 
time and M{t) are given by Eqs. (|2.44|) and (|2.46|) respectively. 

The main conclusion of this Sec. IIIDI lies in the fact that for the moderate and strong forces the translocation 
exponent changes from a — 2v for moderately long chains to a = 1 + for very large values of N. This is in good 
correspondence with the results of the Monte Carlo investigation [3l| . 



III. NUMERICAL VERIFICATION OF THEORETICAL PREDICTIONS 



A. Model 



In order to verify the predictions of Section |TT] we performed numerical simulations. The model assumes Langevin 
dynamics of a polymer chain which consists of N beads that thread through an octagonal pore through a closely- 
packed wall (membrane). The interaction between the monomers of the chain is modeled by a Finitely Extensible 
Nonlinear Elastic (FENE) springs corresponding to a pair potential 

UpENEinj) = In - ^) ' (3.1) 

where rij is the bond length between two beads and Ro — 1.5 is the maximal bond length. All beads experience 
excluded volume interactions which are modeled by the repulsive part of the shifted Lennard-Jones potential, also 




FIG. 3: (Color online) A snapshot configuration during translocation of the polymer chain. The nanopore is created by 
removing 8 beads, so that an octagonal hole results. 



known as the Weeks-Chandler- Anderson (WCA) potential. This potential Uwca is defined by 



UwcAiri 



4e 



/ \12 













e(r, - r), 



(3.2) 



where Q{x) is the Heavyside-function, i.e., we use a cut-off Vc = 2^^/^cr, implying Uwca — for rij > r^. The 
monomers residing inside the pore experience a constant external force / in the direction perpendicular to the mem- 
brane, which we designate by x. The external force can be implemented by adding a linear potential Uext, whose 
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value is outside the pore and fx, if x is inside the pore region. Thus, / pulls the chain towards the region of positive 
X which we refer to as the trans-side. The equation of motion for the beads of the chain reads 

m-^ ^ -V [Ufene + Ulj + Ue^t) ~ 7^ + '^^it), (3.3) 

where TZi(t) stands for a Brownian random force whose moments obey (TZi {t)) = and {TZ" {ti)TZj (t2)) — 2kBT^S{ti — 
t2)5a,p5i,j The parameter values were set to e = 1.2, a — 1.0, k = 60.0, 7 = 0.73 and were kept fixed during the 
simulations. The temperature had a constant value given by T = 1.2e/kB- These parameters correspond to the model 
used in the ref.f27|. 

The membrane is modeled by a plane of beads whose positions are kept fixed. Eight neighboring monomers are 
removed to obtain an octagonal pore. The interaction between the beads of the chain and the plane is mediated 
through the repulsive WCA potential. All simulations were performed starting from a configuration in which initially 
all beads except the first one were on the cis-side. The chain is fully translocated when all beads have made their way 
to the trans-side. To prevent the chain from retracting to the cis-side, reflecting boundary conditions were imposed 
on the first monomer on the chain by making its size larger than the pore diameter. A typical configuration for a 
polymer with N = 100 in 3D is shown in Fig. [S] 



B. Results 



To check the validity of the blob picture, derived in Sec. |TT1 we start by preparing contour plots of the monomer 
density. We assume that the chain propagates in the positive a:-direction, that is, perpendicular to the membrane; the 
z-direction is parallel to the membrane and equivalent to the y. The contour plots of the density distribution Figs.S^ 
- 14}1 describe the change in the average polymer conformation at different stages of the translocation process. 

In order to gain more insight into the translocation process we have also plotted similar data in a different form. 
Figure [S] demonsrates the density profiles (normalized to unity) in x-direction for different time moments. This 
representation clearly shows that the tensile force transmission (or front propagation) in x-direction (as it was described 
in Sec. IIC) is really visible. Namely, the moving domain grows at the expense of the rear part of the chain which 
is still at rest. After approximate time t = 70000, the non-translocated part of the chain moves as a whole (cf. Sec. 
IICl, where this second stage of translocation was qualified as the "suction into the pore"). 

Figure [5] also clearly shows the presence of "crowding" effect, i.e., a fairly strong compression of the coil on the 
<rans-side in x-direction (i.e. at x > 0). This effect results from the fact that the characteristic time t of a forced 
translocation (which ranges, as we have shown in Sec. Ill C| between n ~ N^'^'^ and T2 ~ N"^^), is by all means much 
smaller than the characteristic Rouse time tr ^ N"^^^^ , that is, the threaded beads on the trans-side fail to disperse 
and form an equilibrated coil conformation in an interval t < th. This crowding effect has been discussed recently 
in more detail [32j and it was shown that immediately after translocation an effective Flory exponent Veff = 0.45 is 
observed, i.e., Veff is smaller than v = 0.59 for a chain in equilibrium. 



1. Scaling of the translocation time 

In Sec. Ill CI we have shown that the time for a driven translocation scales as (r) ~ N°'f^^. For large forces (the 
corresponding condition is given by Eq. (j2.27p ) the translocation exponent a ^ 1 + v k, 1.59 and 5 = 1 (for the 
so-called free-draining or Rouse dynamics). Thus one goes back to the scaling relationship, predicted first by Kantor 
and Kardar 12 1. Very recently [0] the same scaling has been obtained on a basis of the so-called "iso-flux trumpet" 



consideration. In fact, in this case the translocation time is determined by the propagation of tensile force with 
characteristic time ti given by Eq. (I2.23p . At relatively small forces, in contrast, the exponents are predicted to be 
a = 2v 1.18, and 6 = l/v 1.7 ("trumpet" regime), or (5 = 1 ( "stem- trumpet" regime). At small forces the process 
is dominated by stationary suction of the rest of the chain after the force transmission stage. 

The results for the (r) vs. N scaling relationship at different driving forces, derived from our MD-simulation, are 
shown in Fig. |6^. We performed numerical simulations for chains with lengths N = 40,70,100,200,300,500 and 
forces 1, 2, 5, 10, 20. For every set of parameters at least 1000 runs were performed, and most averages were obtained 
from 3000 separate runs. The N = 500 runs took a long time, mostly because relaxation to equilibrium (as a starting 
condition) was very time consuming. The / = 1,2 results for the N = 500 chain are therefore typically averaged 
over 1000 runs. It can be seen from Fig. that the theoretically predicted tendency is correct: the translocation 
exponent a grows with increasing force. Nevertheless, quantitatively the MD-values are systematically smaller than 
the theoretically expected ones: a k, 1.33 for the large force and a ~ 1.06 for the small force. This behavior closely 
corresponds to the results found by the Langevin MD method [25l - [27j . 
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Fig. [6]d shows how the translocation time changes with force. For the long chains the exponent S is very close to 
unity in a good agreement with the theoretical prediction. For a shorter chain 5 ~ 1.37 which is much smaller than 
the possible theoretical prediction within the trumpet scenario: d ss 1.7. 



2. Evolution of the translocation coordinate 



Much information about the translocation process can further be obtained by examining the average value of the 
translocation coordinate {s{t)). The results for (s(t)) at different forces and chain lengths are shown in Fig. [T) It 
can be seen that for a fixed force the behavior is largely independent of the chain length, that is, the exponent /? in 
{s{t)) ~ does not depend on iV. Solely the plateau height which naturally (as for any Brownian motion on a finite 
interval) marks the long time limit {s(t — > oo)) linearly depends on the chain length TV. 

One can also verify from Fig. [7] that the slope characterized by the exponent (3 becomes smaller with increasing 
strength of the driving force /. Namely, f3 « 1.06 for small forces and (3 w 0.85 for larger forces. These results could be 
compared with the corresponding theoretical predictions given by Eq. (|2.34l) and Eq. (I2.30|) (see also Fig. [2t )• This 
basically reflects the same tendency which has already been seen in the (r) oc N" scaling relationship. Apparently, 
the translocation time is defined as (s((t))) = N, i.e., (t)^ ~ N and we have (3 w 1/a. This relationships holds well 
in our MD-simulation. Therefore, the fact that a grows with the force, which we have seen before, correlates with the 
decrease of the exponent (3. 

From Fig. [7] it also easily inferred that during the initial time period {s{t)) remains to a very good approximation 
almost constant. We identify this period with the blob initiation time (see Sec. Ill Bp . i.e., the time which is necessary 
to generate a first blob. In Fig. [7] this characteristic time Tinit manifests itself as a first crossover from a constant value 
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FIG. 5: (Color online) Density profiles (normalized to unity ) as function of s-coordinate. In the time interval < t < 70000 
the front moves toward the left, i.e., the moving domain grows at the expense of the rear part of the chain. At t > 70000 the 
second stage of the translocation, "suction in the pore", sets in. On the trans-side (i.e. a; > 0) the "crowding" effect can be 
seen. 




N f 

FIG. 6: (Color online) (a) The t ys N for different forces /. The figure shows that r ~ A'''*, where a depends on the force. 
For small forces (/ — 1, 2), the value of a is about 1.1, for larger forces (/ = 10, 20), a crossover appears at A''* = N 100. 
For chain lengths N < N* one has a « 1.16 which is close to 2v, whereas for TV > A*'*, a = 1.33. (b) r vs / for chain lengths 
A'' = 40, 100, 200, 300, 500. One finds r oc with 5 ~ 1. For short chains lengths the exponent 6 is somewhat larger than 
1. 



(vifithm 10%) to the scalmg law {s{t)) ~ . The time Tinit is chain length independent but inversely proportional 
to the force / (cf. Eq. (|2.2p ). This conclusion is supported by Fig. [S] where the first crossover time Tinit is plotted 
against the force. 



IV. CONCLUSIONS 



In this work we derive scaling laws for the mean translocation time of a driven polymer chain through a narrow pore, 
based on the ideas of tensile force propagation along the polymer backbone p^.ll6||. Our findings can be summarized 
as follows: 
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FIG. 7: (Color online) The average translocation coordinate (s) as function of time for 4 different chains lengths A'^ = 
100, 200, 300, 500 and 5 different forces F — 1,2, 5, 10, 20. The force increases from right to left and the chain length increases 
from bottom to top. The translocation exponent /3 decreases for strong forces. 
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FIG. 8: (Color online) The time needed to create the first blob rinit as a function of force for three dilferent chains lengths. It 
is clear that this time is independent of the chain length and decreases approximately as rinit ~ 1//. For large forces the force 
dependence appears to level off. 



1. The translocation starts with the formation of initial Pincus blob (i.e., the first blob is generated immediately 
at the pore opening) . The characteristic time of the blob initiation is given by Eq. (|2.2p . Our MD-simulation 
results essentially support the scaling prediction Tinit ^ 1/ f. 

2. The initiation is followed by a tensile force transmission along the chain backbone which is governed by the local 
balance of driving and drag forces. For forces in the interval ^ af /ksT < 1 this leads to a trumpet regime 
(see Fig. [1]) . The corresponding translocation time is given by Eq. (j2.26p where the first term, which corresponds 
to the force transmission characteristic time, prevails under the condition Eq. (|2.27p . As a result, depending on 
the force strength and on chain length (for Rouse dynamics), one expects a crossover from 'M'^^ j j'^l^ to N^^^ / f, 
i.e., the translocation exponent a grows with increasing force / from a « 1.18 to a k, 1.59. 

One should note that the MD-simulation findings yield systematically smaller values for the translocation ex- 
ponent: a « 1.06 and a ~ 1.33 for weak and strong forces res pec tively. This is in agreement with other 
simulation results [25l - [27| and differs from findings reported in ref. |28j . The mentioned slight overestimation of 
the translocation exponent a by the theory may be due to the role of fluctuations which are not accounted for 
within the quasi-static approximation used in this work. As for the force scaling, the theory predicts that the 



16 



scaling law ranges between 1/ f^-^^ and 1/f with the chain length increasing. Again the MD-simulation gives 
smaller exponents: l/f^-^'^ and 1/ f^'^^ respectively. 

3. Under the conditions, given by Eq. p.27p . the number of translocated segments M(t) changes as M{t) — 
coilatf'^'^'"'^ (see Eq. i.e., in the scaling law M{t) - {JtY the exponent /3 w 0.63 for relatively large 
forces. Our MD-simulation result gives a slightly larger value fi ~ 0.85. The reason for that is the same as in 
the case of the scaling law (r) ^ A^": recall that p ~ 1/a. 

4. For strong forces the tensile force transmission follows either the "stem-trumpet" (1 ^ af /ksT <^ N^) or the 
"stem" {N^ <^ af /ksT) scenarios (see Fig. [TJd, c). In both cases the translocation time can be estimated as 
(r) = Ci tqN'^+'' /fa + C2 To N^^/fa whcrc the first term dominates under condition N^^" > C2/C1. In other 
words, the translocation scaling exponent grows from a = 2;^ to a = l+v as the chain length N increases. This is 
in agreement with the results of LD-simulation 23] and also with the findings of Monte Carlo investigations 31]. 

We believe that the present approach leads to a consequent and physically plausible theory of driven translocation 
dynamics even though the theoretical predictions regarding scaling coefficients do not fully agree with data from 
computer experiments. The theory is based on the assumption that drag and driving forces equal one another during 
the process of chain threading through the pore. Since at the present level of theoretical treatment fluctuations are not 
taken into account, further developments of the theory should try to incorporate fluctuations too. Also the analytical 
model totally ignores additional resistance forces which are imposed by the pore and by the over-crowded beads at 
the pore exit (as it can be seen from Fig. . These forces should be included in the general force balance given by 
Eq. (j2.5p . On the other hand, these resistance forces make the translocation dynamics even more sluggish, so that one 
can not explain thus the observed decrease of the translocation exponents as compared to the theoretical prediction. 
In principle, one could extend the present consideration on the case of attractive pores where entering monomers may 
be captured and reside for a while in the pore. Such studies have been carried out recently [ss!, 'si']. One should note, 
however, that this effect does not alter the scaling laws for translocation even though the some additional details 
of the process become thereby important. Evidently, more work is therefore needed until a full understanding and 
entirely satisfactory description of translocation dynamics is achieved. 
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Appendix A: Integral and differential forms of the material conservation law 

Here we demonstrate why the calculation, performed in ref. jl6l] . leads to a scaling result for ri which deviates from 
the one given by Eq. (|2.23p in Section III CI In contrast to the integral form of the material balance equation, given 
by Eq. (|2.14l) and used in the present work, in ref. [H, HBl one uses its differential form (cf. Eq. (5) in ref. (iq): 

As before, p{x) ~ g{x)/^{x)^ is the local density and S(a;) ~ C(^)^ stands for the cross-sectional area so that 
x{x) = p{x)Y^{x) can be seen as the linear density. 

One may readily see that the differential representation of the material conservation, Eg. dAll) . suffers from the 
problem that the linear density x(a;)|^^„x(t) = [^(^))/'^]^^''^^/"L=-x(t) '^^^^ exist because it diverges! Indeed, 
it is evident from Eq. (|2.12p that £,{x) at x = —^{t) has an integrable singularity. The physical reason for such 
behavior is clear: the force at the free boundary is zero, therefore, the tensile blob size goes to infinity (cf. Eq. (|2.10l) . 
Because of this, one should use from the very beginning the integral form Eq. (I2.14[) and may turn to the differential 
equation for X{t) (cf. Eq. (|2.20p ) only after calculation of the involved integral. 

Instead, Sakaue [H, Eal introduces a "cut-off" x = —X{t) + ^x{t), where £,x{t) is the "largest blob" given by 
Eq. (|2.13p . The implementation of the "largest blob" £,x{t) as well as the relationships given by Eq. (|2.1ip (at 
X{t) :s> ksT/ f) and Eq. (|2.15l) in Eq. ljAip leads eventually to a differential equation for the front location X{t) (in 
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dimensionless units): 



To 



UJ-\-l £UJ-\-2~-Z 



2-2 



dX{t) 
dt 



(A2) 



where X{t) = X{t)/a, fa = af/kBT, a; = (1 — — 2)/v{z ~ 1), and cq is a numerical coefficient. The solution of 
this equation reads (of. Eq. (|2.22l) ') 



t = t^+ Tof:+'-^xitr+' \i - c^/x{trra 

Finally, by equating X{ti) = N, Sakaue obtains the characteristic time 



Tl - To 



jU+2 — z jiq(u+2)ii 



(A3) 



(A4) 



which is differs from our result, given by Eq. (|2.23l) . In the case of Rouse dynamics (z = 2 + 1/;^), the scaling law 
Eq. (1X4)) leads to a = (1 + + 2v'^)/{l + v) k. 1.43 and 5 + 1.26. Despite the fact that this value for 

the a-exponent correlates better with our MD-findings, one may doubt if the calculation based on the "cut-off" trick 
is trustworthy. Actually, this "cut-off" trick makes the "trumpet" effectively shorter so that Sakaue's characteristic 

time Tl - is smaller than ours n ~ N^-^^/f. 

Contrary to Sakaue's calculation, we use the integral form of the material conservation law Eq. (|2.14p where the 
integral term can be calculated exactly (without resorting to any "cut-off" trick). By making use of the closure 
relation, Eq. (|2.15p . and the relationship for the translocated monomers, M{t) Eq. (|2.17p , we derive a differential 
equation for X{t), Eq. (|2.20l) . which differs from the corresponding Eq. (jA2l) . 
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